Topological Transition in a Non-Hermitian Quantum Walk 
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We analyze a quantum walk on a bipartite one-dimensional lattice, in which the particle can 
decay whenever it visits one of the two sublattices. The corresponding non-Hermitian tight-binding 
problem with a complex potential for the decaying sites exhibits two different phases, distinguished 
by a winding number defined in terms of the Bloch eigenstates in the Brillouin zone. We find that the 
mean displacement of a particle initially localized on one of the non-decaying sites can be expressed 
in terms of the winding number, and is therefore quantized as an integer, changing from zero to one 
at the critical point. This problem can serve as a simplified model for nuclear spin pumping in the 
spin-blockaded electron transport regime of quantum dots in the presence of competing hyperfine 
and spin-orbital interactions. The predicted transition from pumping to non-pumping is topological 
in nature, and is hence robust against certain types of noise and decoherence. 



A quantum system is said to exhibit a topological tran- 
sition when it features several phases, characterized by 
a topological invariant that takes on different quantized 
values in each of these phases [l[ . A classic example of a 
topological transition occurs in the quantized Hall effect, 
which can be linked to the Chern invariant |^ , defined in 
terms of the system's single-particle wavefunctions with 
quasi-periodic boundary conditions. Because the Hall 
conductance is proportional to the Chern invariant, and 
because quantization of the latter is of a topological na- 
ture, the quantized Hall effect is universal across sam- 
ples of varying size, shape, or composition, and is robust 
against many types of disorder. Another example of a 
similar nature is encountered in adiabatic transport 

Here we present a model exhibiting a new type of 
topological transition in a system described by a non- 
Hermitian Hamiltonian. We consider a quantum walk on 
a bipartite one-dimensional (ID) lattice, from which the 
"walker" (particle) can decay whenever it resides on the 
sites of one of the sublattices (see Fig[T^) . Due to hop- 
ping between sites, a particle initially localized on any of 
the non-decaying sites at time t = Q will eventually decay 
from the system as i — > cxd. 

Surprisingly, we find that the average displacement 
of the particle during the course of its decay. Am = 
"^Prn, IS exactly quantized as an integer (0 or 1 unit 
cells), where Pm is the probability distribution for de- 
cay from different sites (see FigHJ;). As in the case of 
the quantum Hall conductance, this quantization results 
from an underlying topological structure; in this case it 
is the winding number of the relative phase between two 
components of the Bloch wave function. Using the topo- 
logical origin of this phenomenon, we are able to show 
that the quantization is insensitive to parameters and is 
robust against certain types of noise and decoherence. 

In recent years, non-Hermitian quantum mechanics in 
one dimension has found applications to a variety of 
different problems, such as vortex matter quantum 
chaosQ, random lasersjl], population biology[3], and 
others (see Ref.0 and references therein). Much of the 
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FIG. 1: Setup of the model, a) Each unit cell m contains two 
sites L (open circles) and R (filled circles), with each 7i-site 
connected to an external decay channel. Intracell (wavy lines) 
and intercell (straight lines) tunneling occur with amplitudes 
V and respectively, b) Energies of the L and 7? sites. Due to 
decay, the ii-site energy obtains an imaginary part, er = er — 
1^,7/2. c) Schematic distribution of local decay probabilities 
{Pm} used to calculate the displacement (|3J. 

interest in these ID problems was triggered by the idea 
that an Anderson localization transition can occur in dis- 
ordered transport with an imaginary vector potential 
In contrast, our problem is translationally invariant; the 
transition results from competition between two pro- 
cesses, intracell and intercell hopping, which occur with 
amplitudes v and v' (see FigUJi). 

In this paper, our motivation is to provide a simple 
model of nuclear spin pumping in spin-blockaded double 



quantum dots 0, [iC 



ll| in the presence of competing 
effects of the hyperfine and spin-orbital interactions, as 
in Ref.[ll[. In the DC transport regime, an electron, first 
loaded into a triplet spin state \L), makes a transition 
to a singlet spin state \R), which is broadened due to its 
couphng to the drain lead [13] (see Fig[T)3). Such transi- 
tions require an electron spin-flip, which can be mediated 
by either the hyperfine interaction or the spin-orbital in- 
teraction, denoted in Fig[T]by the amplitudes v' and v, 
respectively. The hyperfine process is accompanied by a 
change of the z-projection of nuclear spin. Am = ±1 for 
an L-state of the type T±, whereas for the spin-orbital 
process Am = 0. Without loss of generality, here we 
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focus on the case of transport through a T+ state, pic- 
tured in Fig[T] The resulting coherent dynamics in the 
combined Hilbert space of electron and nuclear degrees 
of freedom is thus described by a quantum walk like that 
shown in Fig[T^. The topological transition, which is ac- 
companied by the formation of a non-decaying dark state, 
leads to a prediction of threshold-like pumping of nuclear 
polarization, along with strong suppression of current due 
to the divergence of dwell time at the threshold. 

Because this problem exhibits a topological phase tran- 
sition, here our main aim is to discuss this interesting be- 
havior from a general point of view; the detailed analysis 
of implications for spin-blockade is postponed to a sepa- 
rate publication |13l] . Here we just note that an ensemble 
of nuclear spins can be described by a single variable m 
in the "gia nt sp in" approximation of constant hyperfine 
interaction IJ], which ignores decoherence arising from 
the more realistic position-dependent interaction [l5|. 

The configuration space of the problem is defined by 
the electron states | L ) , \ R) and by the total nuclear 
polarization, taking integer values — oo < m < cxd. 
Thus the Hilbert space has a tensor product structure: 
span{| m) ^ \ L/R)}. In this basis, illustrated in FiglT^, 
the state of the system | -0 ) is described by the ampli- 
tudes = ( m L I -0 ) and = {mR\ip), and evolves 
according to the equations of motion 
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Without loss of generality, we choose v > Q and v' > 0. 
The on-site energy in = er — ih-y/2 for the i?-states has 
an imaginary part that accounts for the decay of these 
states with rate 7, while the on-site energy is real. 
Now, suppose the system is initialized to the L-state 
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at time t = 0, and allowed to evolve freely under the equa- 
tions of motion ([T]). Because of translational invariance, 
we can equivalently start anywhere on the L-sublattice. 
Under the dynamics ((T]), the wavepacket describing the 
quantum walker spreads throughout the lattice and leaks 
out through its components on the i?-sites, decaying com- 
pletely as t 00. What is the average displacement 
achieved by the particle before leaking out? More pre- 
cisely, given the ability to detect the site m from which 
the decays occurs, and thereby measure the decay prob- 
ability distribution Pm (see FigHJ:), we would like to find 
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Although (Am) can be obtained from an explicit calcula- 
tion involving the system's time evolution operator, here 
we will pursue a less direct but more rewarding approach 
that helps to uncover the topological structure behind 



the solution. The result is supported by numerical sim- 
ulations, which also allow us to test various features of 
the model such as its robustness against decoherence. 

As a first step in the calculation of (Am) , we note that 
the norm of a quantum state | if) ) evolves according to 
-^(iplifj) ^ iitp \ {H^ — H) \ ^ )■ For Hermitian systems, 
W = H and ^( | ) =0. However, our system is non- 
Hermitian due to the complex energy in, and, as seen 
from the equations of motion ([T]), decays according to 
^(V'lV') = — SmTlV'mP- The decay is thus described 
as a sum over local terms accounting for the decay from 
each site of the lattice, Eq.dS]), with Pm = 1- 

It is beneficial to pass to the momentum representa- 
tion, -0^ = § dk e^^"^ ijj^ , where the integral is taken 
over the Brillouin zone — tt < k < n. Due to the trans- 
lational invariance of the system ([1]), the equations of 
motion in the Fourier representation break up into 2x2 
blocks, one for each momentum k: 



dt 



0f 



(4) 



with Afc — V -\- w'e . The two-component wave func- 
tion for each k evolves independently of the others, with 
Pk{t) = |'0fc(i)P + iV'fc the probability density to 

find the system with momentum k at time t, decaying as 

dtPk = -^\i^m\'- ^ 

Writing m as a derivative with respect to k via m tp^ = 
~TiT § '^^ Ik (e*''™) "^k ^"^^ integrating by parts to move 
the derivative onto ij]^, we bring Eq.Q to the form 
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Next, we use the polar decomposition i{j^{t) — 
Uk{t)e''^>'^*\ where Uk = \'<P§{t)\ and 0k = arg{'0^(t)}. 
We assume that Uk{t) > for all t > 0, which follows 
from Eq.Q after some algebra Using the fact that 
f dkukdkUk = is an integral of a total derivative over 
a closed contour, we rewrite Eq.([5]) as 



(Am) = -1 dt j '^\uk{t)\ 

dpk d0k 
dt dk ' 
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dk 
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(7) 



where we replaced — 7|ufe(t)|^ by dtPk in Eq.®. With 
the help of integration by parts in the integral over i, the 
time derivative can be moved from pk onto dkOk-, giving 
(Am) =Io~ dt§^pk dtidkOk), with 

d0k 



lo = 



dk 
2^ 



Pk 



dk 



(8) 



We will now show that the boundary term Xq provides the 
only non-zero contribution to the integral ([7]). First, we 
use integration by parts on the integral over k to obtain 



(9) 
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As demonstrated below, this integral vanishes because pk 
and dtOk are both even functions of k. 

In order to see that pk and dt9k are even functions, it 
is helpful to view the evolution (jH) within each 2 x 2 fc- 
subspace as the precession of a decaying pseudospin in a 
(complex) magnetic field with z-component £l — £r and 
transverse component of magnitude 2\Ak\ = 2\v v' e^^]. 
Because \Ak\ = l^-fc|, a static rotation about the z-axis 
maps H-k into Hk, with Hk the 2x2 matrix in Eq.(l4]): 



Hk, ipk = arg{Afc}. (10) 



Given that the initial state ^ is oriented along the z- 
axis in pseudospin-space for all fc, in the rotated frame 
pO)) the pseudospin associated with the momentum — fc 
performs the identical evolution to that of the pseudospin 
associated with momentum k. Because the state ([2]) has 
equal magnitude in all momentum sectors, the moduli of 
the k and —k pseudospins are equal for all times, Pk{t) = 
P-k{i)- Furthermore, from (jlOp their phase difference 
is time-independent, 0_fc — Ok — 2ipk, which proves the 
claim. 

To evaluate To, we use the facts that all fc-states are 
initially occupied with equal probability pfc(t = 0) = 1, 
and that the state decays completely, pk{t oo) = 0. 
Substituting these values in Eq.®, we find 



(Am) 



dk del 

2^ 'dk' 



lim 9k{t). 



(11) 



Although Ufe(O) = 0, the limit t ^ 0+ is well-defined. 

Expression (jlip is a surprising result: the expected dis- 
placement of the particle as it spreads out and decays is 
equal to the winding number of the relative phase between 
components of the Bloch wave function. In particular, 
this means that (Am) can only take on integer values. 
Using ijjkidt) = -iAldt/h, we have 6*^ = arg{-iA^}. 
It is thus immediately clear that there are two possible 
situations depending on whether or not Ak = v + v'e^^ 
wraps the origin as k is taken around the Brillouin zone: 
(Am) = 1 (0) when v' > v {v > v'). 

It is perhaps not entirely obvious from this discussion 
that the transition at u = w' is a characteristic of the 
Hamiltonian rather than of the initial state. To clarify 
this point, we examine the eigenstates of Hk and plot the 
ratio of their components ^k = i'k I'^k complex 
plane. As shown in Fig|2^, the winding number about 
the origin changes from 1 at u' > t; to at u' < t;. 

Furthermore, one of the eigenvalues of Hk becomes real 
at the transition v = v' , because Ak vanishes for k — 
TT. This indicates the formation of a non-decaying dark 
state with i^j}^^ — 0; under these conditions, the k = 
TT component of the initial state ^ remains stuck on 
the I L ) sublattice for all t. Dark states formed in the 
nuclear subspace, as found here, can arise as fixed points 
of cooling processes 14 1. 
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FIG. 2: a) The winding number change for an eigenstate of 
the Hamiltonian Q. Here we plot the component ratio = 
ipk /''Pk momentum k for — tt < k < tt. b) The expected 
displacement, Eq.(|3}, after full decay of the initially localized 
state. The quantization of (Am) is topological in nature, and 
is linked to the winding of Ak 



V + w'e'* around the origin. 



As pointed out in Ref.[16[ (see also [17[), dark states 
in quantum dots can result in current suppression due to 
the Dicke effect. In our case, the average decay time 



t-{^\^)dt 



/dk f"" 
— pk{t)dt (12) 



may become very long near the transition (here we used 
A(^2p\ip) = § dk dtPk and integrated by parts). Close 
to the transition v = v' , when |Afe~7r| <C fvy, the dynam- 
ics dl]) yields Pk{t) ~ exp(— Ffei), where Ffe is given by 
Fermi's Golden Rule: Tk = \Ak\''i /[{eL-eB)^ + {h-i /2Y]. 
Substituting these expressions into Ea. (|12p . and using 



the change of variables z = e^^ , we get 
_ (SL - eu? + {h^/2f 



27ri7 



dz 



{vz + v'){v + v'z) 



, (13) 



where the integral is taken over the unit circle \z\ = 1. 
Using the residue theorem, we see that the decay time 
indeed diverges at v' = v as f oc l/\v — v'\. 

Gonspicuously, neither the quantization of (Am), nor 
the discontinuity at t; = w', seem to depend on the values 
of the decay rate 7 or the energies £l/r- Furthermore, 
the analysis leading up to Eq.(fTT|) goes through even if 
7 and £l/_r are made time-dependent. In particular, the 
integral ([9|) still vanishes because the states k and —k 
see identical time dependent effective fields, up to a ro- 
tation pO|) . This suggests, among other things, that the 
sharp transition shown in Figl2]3 survives dephasing due 
to classical noise on the energy levels El and e^. 

To investigate this remarkable indifference to dephas- 
ing, we have performed direct numerical simulations of 
the equations of motion ([T]) up to a fixed time T and 
restricted to a finite chain of 51 unit cells. During each 
time step tn < t < tn + At, we evolve the state forward in 
time and bin the probability of decay from each unit cell. 
In Fig[2^ we show the results for the mean displacement 
(jlip and the decay time f, obtained using the distribution 
Pm (see FiglSjD), and the formula f = J2n IV'(^n)P^*- 

This simulation, showing clear quantization, was then 
altered to investigate the robustness against decoherence. 
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FIG. 3: Results of simulation with finite diain of length 
= 51 unit cells and 7 = el — er = 1- a) Displacement (Am) 
(blue circles) and decay time r (red diamonds). Filled sym- 
bols were obtained by evolving the wave function up to time 
T = 100. Decay slows down near the critical point, where 
longer running time is required (open symbols, T — 500). 
The black dashed line shows 1 — \tp{T)f, used to monitor 
completion of the simulation. Dotted lines are smooth spline 
fits added as an aid to the eye. Incomplete decay due to fi- 
nite running time T and finite length A'' appear as rounding 
of the step. Results of simulation with r2 — 10 for linear 
damping (green boxes) and repeated projective measurement 
(yellow triangles) show that quantization survives L — R de- 
coherence, b) Decay probabilities {P,n} at v/v' = i, 0.85, 9. 
The distribution becomes broad near the transition v = v' , 
while the mean remains quantized. 

We modified the simulation to evolve the system's density 
matrix, adding an exponential damping of the L — R-oS- 
diagonal elements with rate Due to the increased 
time and memory requirements, we could only simulate 
smaller systems over a more sparse sampling of points. 
However, the results do show a relatively well- formed 
step (green boxes), consistent with expectation. Simi- 
larly, the step appears to be robust against a stronger 
form of decoherence where the density matrix is repeat- 
edly projected onto the L and R subspaces at a fixed 
time interval Ar = l/r2 (yellow triangles). 

The transition is not robust against all types of noise, 
however. Any variations of the amplitudes v and v' in 
time will in general broaden and distort the step. Like- 
wise, we do not expect the sharp step to survive pertur- 



bations that break translational symmetry. 

In our motivating example of spin-blockaded transport 
in double quantum dots, the topological transition would 
be manifested as an abrupt change in the nuclear spin 
pumping rate as the relative strengths of spin-orbit and 
hyperfine matrix elements are varied, for example, by 
tuning gate voltages to change the electrostatic poten- 
tial felt by the electron. Because pumping thresholds 
in various parameters are ubiquitous in such systems, 
to unambiguously identify a pumping threshold with the 
topological transition discussed here one must correlate 
the appearance/disappearance of nuclear spin pumping 
with a decrease in the current through the system result- 
ing from the diverging dwell time f, which is not expected 
for pumping thresholds of other origins. 
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